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TECHNICAL  NOTE  2878 

COMBINED  EFFECT  OF  DAMPING  SCREENS  AND  STREAM 

CONVERGENCE  ON  TURBULENCE 
By  Maurice  Tucker 


SUMMARY 

An  analysis  is  presented  of  the  combined  effect  of  a  series  of 
damping  screens  followed  by  an  axisymmetric -stream  convergence  (or  diver¬ 
gence)  upon  the  mean-square  fluctuation-velocity  intensities,  scales, 
correlations,  and  one -dimensional  spectra  of  a  turbulence  field  con- 
vected  by  a  main  stream.  The  treatment  is  restricted  to  negligible  tur¬ 
bulence  decay  and  linearized  by  postulating  small  fluctuation  velocities 
and  velocity  gradients,  and  absence  of  viscosity  except  as  simulated  by 
the  idealized  screen  action.  Compressibility  of  the  main  stream  is 
allowed  for  during  passage  through  the  contracting  section.  The  density 
fluctuations  associated  with  the  turbulence  field  are  regarded  as 
negligible. 

Numerical  results  for  the  statistical  quantities  describing  the 
turbulence  field  downstream  of  a  screen-contraction  configuration  are 
obtained  for  the  case  of  upstream  isotropic  turbulence.  The  action  of 
the  damping  screens  and  the  stream  convergence  is  to  distort  this  ini¬ 
tially  isotropic  field  into  a  field  of  turbulence  symmetric  about  the 
longitudinal  direction  with  the  lateral  fluctuation  velocities  greater 
in  magnitude  than  the  longitudinal  velocities. 

An  approximate  method  of  taking  into  account  the  effects  of  turT 
bulence  decay  upon  the  mean-square  fluctuation  velocities  obtained  for 
the  case  of  negligible  decay  is  presented.  This  method  of  correction 
together  with  the  tabulation  of  fluctuation-velocity  ratios  over  an 
extensive  range  of  conditions  should  prove  useful  for  engineering 
applications . 


INTRODUCTION 


The 

settling 
attain  a 


use  of  fine-mesh  or  damping  screens  located  in  a  low-speed 
chamber  followed  by  a  contracting  passage  (entrance  cone)  to 
low -turbulence  test-section  flow  is  well  known  from  the 
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qualitative  standpoint.  Dryden  and  Schubauer  (reference  l)  have  pre¬ 
sented  experimental  data  regarding  the  combined  effect  of  screens  and  a 
contraction  on  the  intensity  of  turbulence.  Existing  theoretical  studies 
are  confined  to  either  the  effect  of  the  screens  or  of  the  stream  con¬ 
traction  on  turbulence.  Taylor  and  Batchelor  (reference  2)  have  obtained 
the  effect  of  a  damping  screen  located  in  a  constant -area  passage  upon 
a  triple  Fourier  integral  representation  of  a  turbulent  field.  The 
effect  of  a  contraction  upon  a  similar  representation  is  analyzed  in 
reference  3.  In  both  references  2  and  3  initial  isotropy  is  postulated 
in  order  to  obtain  numerical  results. 

The  analyses  of  references  2  and  3  indicate  that  in  the  absence  of 
decay  effects  (dissipation  and  mixing)  an  initially  isotropic  turbulence 
field  -will  be  distorted  into  a  field  of  turbulence  axisymmetric  about 
the  mean  flow  direction  upon  passage  through  either  a  damping  screen  or 
an  axisymmetric  contraction  (contraction  with  all  cross  sections  simi¬ 
lar)  .  An  analysis  of  axisymmetric  turbulence  is  given  in  reference  4. 

In  conventional  wind-tunnel  configurations,,  turbulence  that  is  initially 
isotropic  will  thus  have  been  distorted  into  axisymmetric  turbulence 
after  passage  through  the  first  of  the  several  damping  screens  and  will 
remain  axisymmetric  while  traversing  the  remaining  screens  and  the  fol¬ 
lowing  contraction.  Inasmuch  as  the  expressions  obtained  in  reference  3 
for  the  downstream  mean-square  velocity  fluctuations  require  that  the 
turbulence  upstream  of  the  contraction  be  isotropic,  the  results  of  ref¬ 
erences  2  and  3  cannot  be  combined  in  any  simple  manner  to  obtain  the 
joint  effect  of  screens  and  a  contraction  on  turbulence  that  is  ini¬ 
tially  isotropic. 

* 

t 

The  present  analysis  treats  the  combined  effect  of  a  series  of  N 
(symbols  are  defined  in  appendix  A)  identical  damping  screens  and  a 
downstream  axisymmetric  contraction  upon  the  longitudinal  and  lateral 
turbulence  velocity  fluctuations,  scales,  correlations,  and  spectra  of 
a  turbulence  field  described  by  a  triple  Fourier  integral.  The  config¬ 
uration  is  shown  schematically  in  figure  1.  Although  compressibility  . 
of  the  main  stream  is  allowed  for  during  passage  through  the  contraction, 
the  density  fluctuations  associated  with  the  turbulence  are  regarded  as 
negligible.  The  assumption  of  small  turbulent  velocity  fluctuations  and 
velocity  gradients  together  with  the  postulated  absence  of  viscosity,  as 
in  references  2  and  3,  implies  the  absence  of  turbulent  decay  processes 
and  linearizes  the  governing  equations  for  both  the  screen  and  con¬ 
traction  effects. 

After  a  discussion  of  the  spectrum  concepts  used  in  the  present 
analysis,  the  preliminary  portions  of  the  analysis  which  borrow  from  the 
results  of  references  2  and  3  are  concerned  with  the  effect  of  a  screen 
and  of  a  stream  contraction  upon  a  representative  wave  or  Fourier  com¬ 
ponent.  Briefly,  the  screen  affects  only  the  amplitude  vector  of  the 
wave;  the  contraction  acts  to  change  both  the  amplitude  and  wave-number 
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vectors.  In  view  of  the  linearized  analysis  and  the  resulting  absence 
of  modulation  or  mutual  interference  between  the  array  of  plane  waves 
making  up  a  field  of  turbulence,  the  correlation  tensor  is  developed 
from  the  results  obtained  for  a  typical  wave.  The  spectral  tensor  is 
obtained  as  the  Fourier  transform  of  the  correlation  tensor.  Turbulence 
velocity  and  scale  ratios  obtained  from  the  spectral  densities  (diagonal 
components  of  the  spectral  tensor)  axe  then  given  in  tabular  form  for 
the  condition  of  upstream  isotropic  turbulence.  The  one -dimensional 
spectra  and  the  correlation-coefficient  curves  for  a  special  case  of 
upstream  isotropic  turbulence  are  also  determined.  An  approximation  for 
taking  into  account  decay  effects  is  suggested.  This  investigation  was 
conducted  at  the  NACA  Lewis  laboratory. 


ANALYSIS  FOR  NEGLIGIBLE  DECAY 

Spectral  Representation  of  Turbulence 

Turbulence  is  often  regarded  as  an  assembly  of  eddies  of  randomly 
varying  size  and  intensity.  The  present  analysis  treats  the  turbulent 
field  as  a  spectrum  of  plane  sinusoidal  waves  with  all  possible  wave¬ 
lengths,  wave-front  orientations,  and  phases.  This  superposition  provides 
the  necessary  three-dimensional  character  to  the  turbulence  representa¬ 
tion.  Large  eddies  thus  are  represented  by  waves  of  large  wavelength 
(small  wave  number).  The  fluctuation-velocity  components  q^(r  =*  1,  2,  3) 

are  represented  at  a  given  instant  by  the  triple  Fourier  integral 

00 

Qy  (k)e  dk1dk2dk3  (l) 

where  x  is  a  position  vector,  Qy  a  wave -amplitude  vector  (reference  3) , 
and  k  a  wave -number  vector  normal  to  the  wave  front.  In  order  that 
the  wave  amplitude  vector  Qy  be  finite,  the  field  of  turbulence 
described  by  equation  (l)  is  assumed  to  occupy  a  bounded  region  and  to 
vanish  everywhere  outside  this  region.  For  the  case  treated  herein  in 
which  the  fluctuation  components  are  related  by  the  incompressible -flow 
form  of  the  continuity  equation 

220yky  =  0  (2) 

T 

the  plane  waves  of  equation  (l)  are  transverse.  In  the  summation  of 
equation  (2)  the  index  y  covers  the  range  of  values  1,  2,  3. 
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In  order  to  obtain  the  spectral  tensor  and,  in  turn,  the  mean- 
square  velocity  fluctuations,  it  will  be  convenient  to  discuss  first  the 
correlation  tensor  and  indicate  its  relation  with  the  spectral  tensor. 
The  correlation  tensor  Rys(£)  is  defined  as  the  spatial  mean  value  of 
the  product  of  the  velocity  component  q^  at  x  and  the  velocity  com¬ 
ponent  qg '  at  £'=£  +  £  as  £  varies  and  the  separation  vector  £ 
of  the  two  points  remains  fixed  during  the  averaging.  If  it  is  assumed 
that  the  field  of  turbulence  is  homogeneous  and  statistically  steady 
and  that  the  field  is  confined  to  a  parallelepiped  of  edges  2D^,  2D2, 


2D3  and  vanishes  everywhere  outside,  the  space  average  is  derived  in 
reference  3  as 

OS 


(£) 


lim 

T  ->oo 


\(^)V(-)e  dkldk2dk3 


—  CD 

where  Qg*(£)  is  the  complex  conjugate  of  Qg(k)  and  T  is  the  volume 

*®1^2®3  P^allslepipe^’  The  expression  lim  — Qy(k)Qg*(k) 

1  is  equivalent  to  the  spectral  tensor  defined  in  reference  5  as 

the  Fourier  transform  of  the  correlation  tensor  Ry-g,(£) 


or 


^8(1) 


00 


—00 


3 

rrs(S  -  6*  (*)  '  (3) 

A  knowledge  of  the  spectral  tensor  permits,  as  will  be  shown,  determi¬ 
nation  of  the  various  statistical  quantities  describing  a  turbulence 
field.  Equation  (3),  which  relates  the  spectral  tensor  to  the  wave- 
amplitude  vector  obtained  for  a  typical  Fourier  component  in  the  absence 
of  any  modulation  effects, is  thus  basic  to  the  present  analysis. 

For  isotropic  homogeneous  turbulence  fields  wherein  the  incompress¬ 
ible  flow  form  of  the  continuity  equation  is  satisfied,  Batchelor  (ref¬ 
erence  5)  has  shown  that  the  spectral  tensor  can  be  written 

’  0(k)  (k\&  -  VO  (4a) 

where  k2  =  k^  +  k22  +  kg2,  =  1  for  y  =  5,  and  =  0  for 

r  4  5. 
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In  matrix  form 


i  2  ,  2 

k2  +  k3 


rr5(^)  -  G(k)  -k-j^kg 


-klk2 
kl2  +  k32 
"*2*3 


-*2*3 


V  +  Y 


It  is  clear  from  the  definition  of  the  correlation  tensor  that  for 
r  £=  0  the  diagonal  elements  of  the  tensor  yield  the  mean-square  velocity 
fluctuations.  In  terms  of  the  corresponding  elements  of  the  spectral 
tensor  (energy  spectral  densities) 


<¥  - 


r^T(k)  dk-jdkgdkj 


The  mean-square  velocity  fluctuations  of  equation  (5)  refer  to 
spatial  averages.  Hot-wire  instrumentation  used  to  obtain  these  fluc¬ 
tuations,  however,  provides  only  time  averages.  Taylor  (reference  6) 
was  able  to  show  that  the  spectrum  of  the  velocity  fluctuations  in  time 
is  the  Fourier  transform  of  the  spatial  correlation  function.  Taylor's 
hypothesis  (reference  7)  that  the  main  stream  carries  along  the  pattern 
of  a  weak  field  of  turbulence  unchanged  past  the  point  of  measurement 
permits  analysis  of  the  hot-wire  output  signal  in  the  form  of  a  one- 
dimensional  spectrum  defined  in  the  equivalent  of  spatial  terms.  The 
relation  between  the  one -dimensional  spectral  densities  F^  and  the 
three-dimensional  spectral  densities  is  easily  shown  by  writing 

equation  (5)  as 


%  = 


Vy(*)  ‘^2^3  ^i  —  /  Fy(k]_)  dk]_ 


The  various  statistical  quantities  which  characterize  a  field  of 
turbulence  may  be  obtained  from  the  one -dimensional  spectral  densities 


as  discussed  in  reference  8.  Noting  that 
the  correlation  coefficients  are  given  by 


*r2Jo 


Fy  ( k]_)  <%  -  1; 


.  .  RYy(rl^°^°)  1 

v  % 


Fy(kl)  cos  k^r-L  dkq 
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The  Fourier  transform  relations  yield 

2Z2  r* 

= — - —  /  Ry  cos  k1r1  dr-|_  (8) 

'-J  0 

Two  sets  of  characteristic  lengths  are  customarily  defined  for  a  tur¬ 
bulence  field.  The  turbulence  microscales  (mean  lengths  weighted 

in  favor  of  the  small  eddies  which  are  responsible  for  the  greater  part 
of  the  viscous  dissipation)  are  given  by 


k12Fr(k1)  dkx 


(9) 


The  turbulence  scales  (mean  lengths  representative  of  the  average 

size  of  all  the  eddies)  are  obtained  as 


(10) 


This  physical  meaning  for  the  scale  of  turbulence  is  only  applicable 
when  Ry  >  0  as  r-]_->oo. 


Plane-Wave  Analysis  for  Damping  Screens 

The  preceding  equations  indicate  that  the  statistical  quantities 
describing  a  field  of  turbulence  may  be  obtained  from  the  spectral 
tensor  of  equation  (3),  which  is  presented  in  terms  of  the  plane-wave 
amplitude  vectors  Qy-(k) .  The  assumptions  of  small  turbulent  velocity 
fluctuations  and  of  inviscid  flow,  with  regard  to  both  the  main  stream 
and  the  turbulence  field  convected  by  the  main  stream,  linearize  the 
equations  which  govern  the  action  of  the  screens  and  of  the  contraction. 
In  the  resulting  absence  of  any  modulation  or  interaction  effects 
between  waves,  the  analysis  is  simplified  by  first  treating  the  effect 
of  a  screen  and  a  stream  convergence  (or  divergence)  upon  a  represent¬ 
ative  plane  wave.  Superposition  is  then  used  to  obtain  the  combination 
of  these  effects  upon  the  complete  assembly  of  plane  waves  which 
describes  the  turbulent  field. 

The  action  of  a  fine-mesh  or  damping  screen  on  a  disturbance  con¬ 
vected  by  a  low-speed  uniform  stream  may  be  characterized  by  two  param¬ 
eters  K  and  a.  The  parameter  K  is  defined  in  terms  of  the  pres¬ 
sure  drop  /IP  required  to  drive  fluid  of  density  p  and  velocity  U 
through  the  screen 
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The  parameter  a  which  takes  into  account  the  side  force  per  unit  area 
was  introduced  by  Taylor  in  reference  9  and  relates  the  angles  of  flow 
incidence  and  flow  emergence  \|f£  shown  in  figure  2.  It  has  been 

shown  experimentally  that  the  ratio  tan  \[rg/tan  tp  tends  to  a  finite 

limit  a  as  i|fp,  which  is  usually  very  small,  tends  toward  zero.  For 
incompressible  flow  the  continuity  equation  requires  that  the  longitu¬ 
dinal  velocity  component  be  unchanged  after  passage  through  the  screen. 
From  kinematical  considerations,  at  the  screen  the  ratio  of  downstream 
to  upstream  lateral  velocity  components  equals  a  for  small  values  of 
the  flow  incidence  angle  \|r-^. 


As  in  reference  2  the  uniform  stream  is  regarded  as  incompressible 
and  inviscid  throughout  the  constant-area  settling  chamber  in  which  the 
screens  are  located  (station  A  to  station  B  of  fig.  l) .  A  screen  will, 
in  general,  decrease  turbulent  motions  of  larger  scale  than  the  mesh 
size  and  introduce  turbulence  of  smaller  scale.  In  the  analysis  the 
damping  screens  are  assumed  not  to  generate  any  wake  turbulence,  which 
implies  that  the  screen  mesh  size  and  wire  diameter  are  very  small  rel¬ 
ative  to  the  scale  of  the  upstream  turbulence.  Far  upstream  of  the 
screen,  at  station  A,  a  single  plane  wave  carried  along  by  the  main 
stream  of  velocity  U  in  the  xp-direction  will  be  designated 


A  0i(k.  x-k-jUt) 


Coordinate  axes  are  fixed,  with  the  origin  located  at  the  screen  and  the 
positive  xp-axis  pointing  downstream.  It  is  shown  in  reference  2  on  the 
basis  of  a  steady-state  disturbance  analysis  that  far  downstream  of  the 
screen,  at  station  B,  the  wave  is  transformed  to 

~  B  m  B  iU'X-kpUt) 

q.r  =  Qy  e 


In  order  to  satisfy  conditions  at  the  screen,  it  is  necessary  to  postu¬ 
late  disturbance  fields  upstream  and  downstream  of  the  screen  which  are 
induced  by  the  screen.  These  disturbance  fields  attenuate,  vanishing 
at  stations  A  and  B.  Taylor  and  Batchelor  represent  these  induced 
velocities  in  terms  of  potential  flows.  With  the  velocity  components 
u,  v,  w  of  figure  2  designating  the  combined  effect  of  the  turbulent 
velocity  fluctuations  and  the  induced  velocities,  the  following  condi¬ 
tions  are  imposed  at  the  screen  (xp  =  0) 


(u) 


B 

Xp=0 


(u) 


A 

Xp=0 
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(v^)x1=0  “ 

The  root-mean-square  fluctuation  velocities  are  taken  to  be  small  rela¬ 
tive  to  the  stream  velocity  so  that  the  equations  of  motion  can  be 
linearized.  A  further  condition  is  imposed  that  the  local  pressure  drop 
across  the  screen  is  determined  by  the  local  longitudinal  velocity  and  the 
screen  pressure-drop  coefficient  K.  The  basic  relations  describing 
this  idealized  action  of  a  damping  screen  on  a  representative  plane  wave 
are  then  obtained  in  reference  2  as 


~  B  ~  A  (0+i)(2a£-iv) 

Qi  ~Qi  (e-i)(2fs+in) 


+ 


iQlAkik2 

tZ 


B.(a-l)2  +  i(v-au) 
L  (p-i)  (20+40 


iQlAkik3 


B  (cl-1)  +  i(p-au) 

.  (P-i)(20+ip) 


(11a) 


( lib) 

(11c) 


where 


£2  s  *22  +  k32, 


2  = 

“g2  ' 


p  =  ( 1+a+K) ,  and  v  s  ( 1+a-aK)  • 


Plane -Wave  Analysis  for  Contraction  Section 

The  main  stream  will  be  regarded  as  compressible  and  inviscid 
throughout  the  contraction  section  (station  B  to  station  C  of  fig.  l) . 

In  the  case  of  supersonic  test-section  flow,  the  term  "contraction"  is 
retained  for  convenience.  As  before,  the  turbulent  field  is  taken  to  be 
incompressible  and  inviscid.  The  contraction  section  has  its  initial 
breadth  and  height  reduced  by  the  factors  Zg  and  23>  respectively, 

while  the  velocity  U(xq)  at  station  B  is  increased  to  ZqU(x]_)  at 

station  C.  A  cubical  fluid  volume  element  of  edge  D  at  station  B  will 
have  been  distorted  into  a  parallelepiped  of  edges  ZpD,  ZgD,  ZjD  upon 

reaching  station  C  (fig.  3) •  The  effect  of  a  contraction  upon  a  tur¬ 
bulent  field  arises  principally  from  changes  in  vorticity  following 
such  distortion  of  the  fluid  elements  passing  through  the  contraction. 

At  station  B  (time  t=0)  in  figure  3,  a  particle  at  distance  x 
from  a  corner  particle  of  a  given  fluid  element  will  at  station  C 
(time  t=t)  be  at  a  distance  X  from  the  corner  particle.  The  coordi¬ 
nate  axes  are  taken  to  move  with  the  main  stream  at  velocity  U(xq) . 

With  the  assumption  of  a  weak  turbulence  field,  the  relative  displace¬ 
ment  of  adjacent  particles  in  a  given  fluid  element  due  to  turbulent 
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mixing  is  taken  to  be  very  much  smaller  than  the  displacement  due  to  the 
contraction.  The  relation  between  x  and  \  is  then  simply 

XT  =  7yXy.  (12) 

I 

With  equation  (12)  the  continuity  equation  for  the  main  stream  in 
Lagrangean  form  provides  the  relation 

alll2l5  =  1  (13) 

where  a  is  the  ratio  of  stream  density  at  station  C  to  stream  density 
at  station  B.  The  product  l represents  the  ratio  of  tunnel  cross- 

sectional  area  at  station  C  to  tunnel  area  at  station  B.  The  parameter 
represents  the  speed  ratio  referred  to  these  stations. 

The  equations  describing  the  changes  in  vorticity  following  dis¬ 
tortion  of  a  fluid  element  are,  from  reference  9: 


c  b  Sxr 

'  -  0  Lj  3xT 


Use  of  equation  (12)  linearizes  these  equations  relating  the  upstream 
and  downstream  vorticities  to 

oyC  =  alyoy®  (14=) 

Upstream  of  the  contraction  at  station  B,  a  single  plane  wave  being 
carried  along  by  the  main  stream  is  designated  at  time  t  =  0  by 

.  =  0yB  (15) 

The  vorticity  at  station  B  is  obtained  from  the  curl  of  equation  (15). 
A  velocity  distribution  at  station  C  compatible  with  equation  (14)  and 
satisfying  continuity,  equation  (2),  is  obtained  in  reference  3  as 


~  r  ~  p  iK-Y 

\  e 


where  the  wave -amplitude  vector  is 
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and  where  the  new  wave -number  vector  K  resulting  from  distortion  of 
the  fluid  volume  element  is  given  by 


ll  l2}  l3 


(17) 


Thus  both  the  wave -number  and  wave -amplitude  vectors  of  a  plane  wave  are 
altered  in  going  through  a  contraction,  whereas  only  the  amplitude  vec¬ 
tor  is  altered  in  traversing  a  screen. 


Equations  (16)  and  (17)  describe  the  effect  of  an  arbitrary  con¬ 
traction  on  a  representative  plane  wave-  For  an  axisymmetric  contrac¬ 
tion  defined  by  the  condition  Zg  =  Z3  (all  cross  sections  are  similar 

but  not  necessarily  circular),  equation  (16)  with  the  aid  of  equation  (2) 
simplifies,  in  expanded  form,  to 


%B  kx2  +  £2 
ll  £kx2  +  £2 


1_ 

l2 


„  B  Qp^kgCl-e) 


Q 


ckp2  +  C2 


QlBkik3d-e) 

€kx2  +  £2 


(18a) 

(18b) 

(18c) 


2  /  2 

where  e  e Zg  /Zp  •  For  an  axisymmetric  contraction,  the  contraction 
parameters  Zp,  Zg,  and  €  may  be  expressed  in  terms  of  the  Mach  num¬ 
bers  at  stations  B  and  C  as  follows : 


'N 


J 


(19) 
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Spectral  Tensors  for  Multiple-Screen-Contraction  Configurations 

Equations  (ll)  and  (18)  describe  the  effect  of  a  screen  and  an 
axisymmetric  contraction  (12=1?),  respectively,  upon  the  amplitude 

vector  of  a  single  plane  wave  typical  of  the  assembly  of  waves 

representing  the  turbulence  field  (equation  (l)).  In  the  Fourier  inte¬ 
gral  q^  corresponds  to  dq^.,  Qy(j&)  to  Cfy  dk-Ldk2'^3?  and  C^Qt)  to 

Qy-  dK^dt^dKj-  Since  at  station  C  (fig.  l)  the  distortion  resulting  from 

the  contraction  transforms  the  wave-number  vector  from  k  to  K_  and 
that  for  axisymmetry  dk-]_dk2dk3  =  Z-j_Z22  (3-ttl<3-*2t3-K3>  equations  (ll)  and 
(18)  yield  "  x  "  °’ 


n  B  __  n  A  (8+i)  (2a|3-iv) 

Qi  -  Qi  (r-T)'C2P+iV) 


+  ^l^A  p(a-l)2  +  i(v-c4x)~l 

t*  L  (P-i)  (26+40  J 


iQlAklk5  |~6(a-l)2  +  i(u-au)~| 
52  L  (6-i)(2p+iq)  J 


,  2_  B 
12  Ql 


+  i* 


.2  »2 

<cki  +  i 


lll2 


lll2 


B  QiBkik2(l-e) 


Q2  + 


+  V 

Bi 


„  B  ,  Qi  kik3(l-e) 

+  2  „2 

ekl  +  C  . 


(20a) 

(20b) 

(20c) 

(20d) 

(20e) 

(20f) 


If  the  fluid  element  volume  T  is  taken  to  be  a  cube  of  edge  D 
at  station  A,  and  hence  at  station  B,  the  volume  will  have  been  dis¬ 
torted  into  a  parallelepiped  of  edges  Z^D,  Z2D,  Z2D  at  station  C  for 

an  axisymmetric  contraction.  The  energy  spectral  densities  which  enter 
directly  into  the  calculation  of  turbulence  fluctuation  velocities  are 
obtained  from  equation  (3)  as 


iim  j^Qr(k)Qr*(k)j  A'B 


(21a) 
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=  Dlima!)  (21b) 

X  2 

With  the  products  QyQy*  from  equations  (20)  formed  and  with  the  use  of 
equations  (21)  and  the  continuity  relations  Qyky  =  0  and  ~ 

the  energy  spectral  densities  may  be  written  as: 


r  B(*) 

ii  - 


4a“k1‘+V2{;‘ 
~~2  272" 

4ki  +p  £ 


r  A(k) 
11  ” 


(22a) 


jj1 22^—^  +  ^33^— ^ 


nB 


-  +  r33h)]A  + 


(V£ 


a 


Vid 


..  2  2,2 
4ki  +  p  £ 


rnA(k) 


(K) 


'k^+C2  X2 


(22b) 

(22c) 


k^d-c)2^  - 


UkxV)2 


( 22d) 


With  the  use  of  equations  (22a)  and  (22c),  the  longitudinal  energy  spec¬ 
tral  density  at  station  C  for  N  screens  in  series  followed  by  an 
axisymmetric  contraction  may  be  expressed  in  terms  of  the  spectral  den¬ 
sity  at  station  A  as 


'1(£L T 

U  W+s2/ 


ta^Wc2^ 
4h12V2d  ) 


(23) 


For  conciseness, equations  (22a),  (22b),  and  (22d)  maybe  written  as 
H-lB  =  AHA,  Vxb  =  +  ZHa,  and  V^c  =  l-j_ V-^6  +  l-p %B,  respectively. 

Then  for  N  screens  in  series,  H^B  =  AH®  ^  =  ANHA,  V^8  =  a2v^_1  + 
and  %°  -  l  The  lateral  energy  spectral  densities  at 

station  C  for  N  screens  in  series  followed  by  an  axisymmetric  con¬ 
traction  may  then  be  grouped  as 


Ixf h  2-a2)j.2)  It-L2 
4k 12+h25Z 


1 


U-eR2 

ek-^+l;2 


4a2k]_2+V2(;2 
4k  -[dn2!;2 


,N-1 


(24) 
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Equations  (23)  and  (24)  relate  the  energy  spectral  densities  downstream 
of  a  multiple -screen-axisymmetric -contraction  configuration  to  the  cor¬ 
responding  upstream  spectral  densities  at  station  A. 


Results  for  Negligible  Decay 

The  solutions  to  be  given  (see  appendixes  B  and  C)  will  now  be 
restricted  to  the  case  of  isotropic  upstream  turbulence.  The  upstream 
energy  spectral  densities  r^A(k)  may  then  be  obtained  from  equa¬ 
tions  (4)  • 


Turbulence  velocity  ratios.  -  As  shown  in  appendix  B,  the  turbulence 
velocity  ratio  or  ratio  of  mean-square  fluctuation  velocities  downstream 
of  a  series  of  N  identical  screens  followed  by  an  axisymmetric  con¬ 
traction  to  the  corresponding  upstream  fluctuation  velocities  is  given 
for  initially  isotropic  turbulence  by 


(^?)A 


'1  JO 


=  (M = a2  ML 
kzf  MA  MA 


5(a2-l)2(V2-a2n2) 


TO  ■*, 

A  sin°  0  d0 

/  2  2  a.2 
(a  -  cos  0) 


_ sin5  8  cos2  9  d6 _ 

(4  cos2S+J  sin2  S)(a2-cos2  6) 2 


,  2—1  ,  .  _4a,2  cos20  +  V2  sin2© 

where  a  =  — —  and  A  = - 2 - 2 - 2~ 

_€  4  cos  0  +  p  sin  0 


A  convenient  approximation  for  equation  (26)  is  presented  later  (see 
equation  (39)). 


For  N  =  1,  equation  (25)  for  the  longitudinal  turbulence  velocity 
ratio  integrates  to 


(j?): 


3a4T}2V2 

4712p2e2(a2-n2)2 


^a2-r|2)(|2  -  a2 )  +  A/i  tanh-1  |)+A2( 


V2-4a2 


where 
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A±  =a2(a2+l)  +  T]2  ( l-3a2)  +  ,+  (a2+l)  (t)2-1)~ 

a2 

A2  =  2(t]2-1)(t12-|2) 


Equation  (26)  for  the  lateral  velocity  ratio  (see  appendix  B) 
integrates  for  N  s  1  to 


For  the  case  of  axisymmetric  contraction  with  the  screen  absent 
(a2=l,  K->  0) ,  equations  (27)  and  (28)  reduce,  respectively,  to 


and 
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which,  in  the  present  notation,  are  identical  with  the  corresponding 
results  of  reference  3.  Similarly,  for  the  case  of  a  screen  and  no 
contraction  (a2->-oo)^  the  results  of  reference  2  are  recovered  in  the 
form 


-  1  + 


3(l-n2)j^l  -  U2-£2)(|  tanh'1 


Punched-card  equipment  was  used  to  obtain  the  turbulence -velocity 
ratios  listed  in  table  I.  For  the  cases  N  =  2,  3,  and  4,  the  integra¬ 
tions  required  for  equations  (25)  and  (26)  were  performed  numerically 
by  use  of  Simpson's  rule  after  changing  the  variable  of  integration 
from  0  to  x  by  applying  the  transformation  x  =  cos  0.  Intervals 
Ax  =  0.01  were  used  in  the  range  0  ^  x  ^0.9;  intervals  Ax  =  0.001 
were  used  in  the  range  0.9|x|1.0.  In~a.ll  computations  the  Mach 
number  Mg  upstream  of  the  contraction  was  taken  equal  to  0.05.  The 
turbulence  velocity  ratios  listed  in  table  I  may  be  corrected  for  values 
of  Mg  other  than  0.05  as  follows:  Values  of  the  parameters  l , 
l22>  and  a^=  — -  for  Mg  =  0.05  and  for  the  desired  value  of  Mg 


are  obtained  from  equations 

,I112)h  ...a  ,  2 


wT 


and  Z2  r== 


(19).  Noting  that  the  quantities 


depend  only  upon  ad  and  K,  the  values 


a2 

l- 


corresponding  to  the  desired  Mg  are 


^  and  known  for  Mg  =  0.05  and 


of  these  quantities  for  the 
obtained  from  table  I.  With 
the  desired  Mg,  the  corrected  velocity  ratios  are  obtained  by  simple 
computation.  The  following  empirical  relations  (reference  l)  were 
utilized  in  obtaining  numerical  results : 


“2  -  (IsJ  for  K  ^ 1 

a2-(pr)  for  K>1 


For  design  purposes,  the  screen  pressure-drop  coefficient  K  may 
be  estimated;  according  to  reference  10,  from  the  solidity  ratio  b, 
where  b  is  the  area  of  the  holes  in  a  unit  area  of  screen,  as 
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For  square -mesh  screen  with  wire  diameter  d  and  mesh  designation  m, 
the  solidity  ratio  as  defined  is 

b  =  (1-md)2 

A  better  agreement  with  the  screen  data  given  in  reference  1  is  obtained 
from 

b  ~  (l-md)7/4 

The  variation  of  the  longitudinal  and  lateral  root-mean-square 
velocity  ratio  with  speed  ratio  Zq  for  a  single  screen  (N=l)  upstream 
of  the  contraction  is  plotted  in  figures  4(a)  and  4(b),  respectively, 
for  selected  values  of  the  screen  pressure-drop  coefficient  K.  The 
results  for  K  =  0,  which  correspond  to  the  case  of  stream  convergence 
or  divergence  in  the  absence  of  any  screen,  are,  of  course,  identical 
with  the  results  of  reference  3.  In  general,  both  the  longitudinal  and 
the  lateral  fluctuation  velocities  downstream  of  the  screen -contraction 
configuration  are  reduced  as  the  screen  parameter  K  is  increased.  The 
somewhat  anomalous  trend  of  the  longitudinal  velocity  ratios  for  values 
of  the  speed  ratio  less  than  2  seems  to  reflect  the  variation  of  the 
auxiliary  screen  parameter  which  approaches  zero  at  K  =  2.76, 

becomes  infinitely  large  in  the  negative  sense  as  K  increases  to  5.28, 
and  becomes  infinitely  large  in  the  positive  sense  as  K  decreases  to 
5.28. 


The  losses  incurred  through  the  use  of  damping  screens  are  propor¬ 
tional  to  the  product  NKU^,  where  N  denotes  the  number  of  identical 
screens  in  series  (multiple  screens)  and  NK  is  the  over -all  screen 
pressure-drop  coefficient.  The  velocity  ratios  for  a  multiple -screen 
arrangement  upstream  of  a  contraction  are  compared  on  the  basis  of  equal 
screen  losses  in  figure  5  for  the  particular  case  NK  =  6-  The  advan¬ 
tages  of  using  a  number  of  screens  in  series  to  attain  a  given  over -all 
coefficient  NK  are  obvious.  An  examination  of  table  I  indicates  that 
the  use  of  multiple  screens  to  attenuate  the  downstream  fluctuation 
velocities  becomes  more  effective  as  the  over-all  coefficient  NK  is 
increased.  The  screen  losses  can  be  reduced  by  decreasing  the  settling- 
chamber  stream  velocity  U^.  Low-turbulence  wind  tunnels  are  generally 
characterized  by  their  many  damping  screens  and  large-cross -sectional- 
area  settling  chambers. 

One-dimensional  spectra.  -  In  accordance  with  equation  (6),  the 
one -dimensional  spectra  at  stations  A  and  B  are  given  by 
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oo 

-JT- 


'^(k)  dk2dk3 


00 


yy  (ii)  d-Kgd-Kg 


As  pointed  out  in  reference  3,  a  comparison  of  the  upstream  and  down¬ 
stream  spectra  on  the  basis  of  the  upstream  longitudinal  wave  number  kq 
is  equivalent  to  a  comparison  of  the  time  spectra  indicated  by  fixed 
hot-wire  probes  located  at  the  corresponding  stations.  Defining  the 


downstream 


spectra  F.^  ( k-|_)  =  Z^F^,2  such 


JFr  (k-^)  dk^  =  /  Fr  ( K-^)  dK-j_,  the  one -dimensional  spectra  at 

0  Jo 

station  C  are  given  by 


[FrC(ki>]N  =  J  4). i 


dkgdk3 


Evaluation  of  equations  (29)  and  (30)  requires  that  thfe  amplitude  func¬ 
tion  G(k)  in  equation  (4)  be  specified.  Compatible  with  the  empirical 
relation  for  isotropic  turbulence  obtained  in  reference  8,  this  function 
is  taken  to  be 


G(k)  = 


(k12+n2+^2)3 


where  the  constants  n  and  H  are  defined  as  n  =  — - -  and 

P  /— \ A  (Lp)A 

H  ,  ■ 


As  shown  in  appendix  C,  the  one -dimensional  spectra  obtained  from 
equations  (4)  and  (31)  maybe  expressed  in  terms  of  a  dimensionless  wave 
number  k-^/n  as  incorporated  in  the  following  parameters: 

s  =  1  +  kq2/ n2 

f  s  s/q2  +  4/p2 

,.o  9  /  p  7  (32) 

g  =  s/|2  +  4a2/v2 

n  2  2 
h  =  - r> - 


a 
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Thus  the  upstream  one -dimensional  spectra  for  this  special  case  of 
isotropic  turbulence  are,  in  dimensionless  form, 


ElAUl/n>  i 

(33) 

FpA(0) 

F2A(kp/n)  _  3s-2 

(34) 

F2A  (0)  s2 

Also,  the  longitudinal  one -dimensional  spectrum  downstream  of  a  single- 
screen-axisymmetric -contraction  configuration  may  be  written  (see 
appendix  C)  in  dimensionless  form  as 


F  p^(  kp/  n) 
FpA(0) 


217 


7  2„2 

ip  p 


„  n  s+h  ,  n 

+  Co  log„  q  +  C 


n  4(s-l) 

3  loSe  — - 

P  s  . 


(35) 


where 


C-,  s 


1_ 

h3 


g(2f -h) 


C2  “ 


P 


gh  h(l+2g)  (l+h)V 

a2(v2-4a2)-  V2 

1 

2fs  f  v2 

_a2(p2-4)  -  p2_ 

J 

J  3sg  g(s+2)  +  2s 

H2( s+h) (g+h) 

p2]\h2  h 

h[a2(p2- 

-4)  -  p: 

_  (s-f)(f-R)(u2-4)2a4 

3  ~  „3r2/  2  .s  212 
f  [a  (p  -4)  -  pj 


The  corresponding  lateral  one -dimensional  spectrum  is 


F2C(ki/n) 

f2a(o) 


2(b-1) 


_  „  „  s+h  „  n  4(s-l) 

E1  +  E2  loge  I"  +  E3  l0ge 


(36) 


where 


* 
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» 


=  a2 (5s -2)  (v2-a2g2)(2s-f)  fg(h-f)  j_  ~] 

1  2s2( s-l)  2„2f2R  a.2fh  L  fh  "  2s  ' 


2|a  f  s 


a  fh 


-Vc 


2a4p2h 


2 ( s+h) 


f(g+h)  -g(h-f)  n2ra2(v2-4a2)-y2]l 

f2h2  v2h2[a2l(i2-4)  -(i2] 


2h(g-f)-fg 


f2h 


E- 


( s+h)[]a2  (y  2  -4a2 


a2h4[a2([i2-4)  -  p2] 


i  -  *5  J  2h  +  55+H  + 


4(a2-l)  (\72-a2p2)h 


a  2  j^a2(|j.2_4)  _  ^2JQ,2(v2_4a2)  _  j 


E*  B 


_  4(s-l)  (v2-a2p2) 


p4f3 


1  + 


L  a2(p2-4)  _ 


The  one -dimensional  spectra  given  "by  equations  (33)  to  (36)  are 
applicable  when  the  amplitude  function  G(k)  has  the  particular  form  of 
equation  (3l) •  Although  these  spectra  are  hot  expected  to  be  valid  for 
the  very  high  wave  numbers  because  of  the  neglect  of  viscosity,  various 
experiments  on  isotropic  turbulence  have  indicated  that  equations  (33) 
and  (34)  provide  a  very  good  approximation  to  that  portion  of  the  actual 
isotropic  spectrum  containing  the  largest  part  of  the  turbulent  energy. 
Equations  (35)  and  (36)  should  furnish  a  similar  approximation  for 
axisymmetric  turbulence.  The  restrictions  given  for  equations  (33)  to 
(36)  do  not  apply  to  the  expressions  for  turbulence  velocity  ratios, 
equations  (25)  and  (26),  for  which  there  is  no  need  to  particularize  the 
spectrum  amplitude  function  G(k) . 


The  downstream  longitudinal  and  lateral  one-dimensional  spectra, 
equations  (35)  and  (36),  are  compared  with  the  corresponding  upstream 
isotropic  spectra,  equations  (33)  and  (34),  in  figures  6(a)  and  6(b), 
respectively,  for  the  following  typical  case:  Mg  =  0.05,  =  2.0, 

K  =  2,  N  =  1.  The  case  K  =  0,  as  obtained  in  reference  3,  has  also 
been  included  for  comparison.  The  scaling  factors  indicated  by  equa¬ 
tions  (B5)  and  (B6)  of  appendix  B  have  been  incorporated  in  the  down¬ 
stream  spectral  ordinates  so  that  the  zero -wave -number  intercept  gives 
the  turbulence  scale  ratio  (appendix  B) . 

The  distortion  in  shape  of  the  longitudinal  spectrum  noted  in  ref¬ 
erence  3  as  a  consequence  of  the  stream  convergence  is  accentuated 
(fig.  6(a))  by  the  presence  of  a  damping  screen  upstream  of  the  contrac¬ 
tion.  This  distortion  is  accompanied  by  a  reduction  in  the  ordinate 

values  by  the  factors  (c^2)  /(,q.q2)  and  (^qq2)  / (lq2)  for  K  =  2 
and  K  =  0,  respectively.  The  downstream  lateral  spectrum  ordinates 
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are  increased  by  the  factors 
B 


Wf/W)‘ 


w;  804 

for  K  as  2  and  K  =  0,  respectively.  The  distortion  in 
shape’ is  relatively  slight  compared  with  the  distortion  noted  for  / the 
longitudinal  spectrum. 


As  may  be  seen  from  equations  (33)  and  (34)  for  the  upstream  iso¬ 
tropic  spectra,  the  longitudinal  and  lateral  spectral  ordinates  have 
maximum  values  at  kp/n  =  0  and  kp/n  =  l//\/3/  respectively.  The 

situation  is  reversed  for  the  downstream  spectra.  Here  the  lateral 
spectral  ordinates  have  maximum  values  at  kp/n  =  0  and  the  longitu¬ 
dinal  spectral  ordinates  at  kp/n  Z  1.4.  Occurrence  of  a  peak  in  the 

spectrum  curve  at  some  wave  number  other  than  zero  is  an  indication 
that  the  correlation  coefficient  may  take  on  negative  values. ^ 


Scale  ratios  and  correlation  coefficients.  -  For  the  scales  of  tur- 
bulence  defined  by  equation  (10),  the  longitudinal  and  lateral  turbulence 
scale  ratios  (ratios  of  downstream  to  corresponding  upstream  scales)  for 
a  screen-contraction  configuration  are  obtained  in  appendix  B  as 


The  scale  ratios  obtained  from  equations  (37)  and  (38)  which  do  not 
require  that  the  amplitude  function  G(k)  of  equation  (3l)  be  specified 
are  listed  in  table  I  for  the  case  of  isotropic  turbulence  at  station  A. 
Typical  results  are  plotted  in  figure  7. 


The  lateral  scale  ratio  (see  fig.  7(a))  approaches  a  constant  value 
of  approximately  4/3  for  values  of  the  speed  ratio  Ip  greater  than  3. 
Measurements  of  the  lateral  correlation  curve  at  a  speed  ratio  near 
unity  which  are  reported  in  reference  11  indicate  that  the  lateral  scale 
is  substantially  unchanged  by  damping  screens.  This  is  in  qualitative 
agreement  with  the  present  result  which  indicates  that  for  7p  slightly 
greater  than  unity  the  downstream  lateral  scale  will  not  exceed  the 


P-1  -kp/n 


"''For  example,  when  F^-(kp)  =  kp  e 
cient  is  obtained  by  using  equation  (7)  as 

Us'  * 


T  (P)  cos  ( '] ?  tan-^-  ar.|J  where 


the  correlation  coeffi- 


F  designates  the  gamma 


function.  For  P  =  1,  R^rp)  is  always  positive;  for  P  >  1,  Ry(rp) 
will  take  on  negative  values  for  particular  values  of  rp. 
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corresponding  upstream  scale  by  more  than  about  20  percent.  Taking  the 
lateral  scale  ratio  equal  to  4/3  leads  to  the  following  convenient 
approximation  for  the  lateral  turbulence  velocity  ratio  from 
equation  (38) : 


(39) 


For  a  given  value  of  the  screen  pressure-drop  coefficient  NK,  the 
longitudinal  scale  ratio  (see  fig*  7(b))  decreases  with  increasing  speed 
ratio  l to  a  minimum  value  at  l ^  =  27*4  (corresponding  to  Mg  —  0.05, 
Mg  -  ^/3)  where  the  contraction  parameter  a2  has  its  minimum  value. 

As  shown  in  table  I,  the  longitudinal  scale  ratio  attains  a  zero  value 
when  the  screen  parameter  u2  =  0  (K  ~  2.76).  This  and  the  occurrence 
of  maximums  in  the  downstream  longitudinal  spectrum  curves  at  nonzero 
wave  numbers  suggest  that  the  downstream  longitudinal  correlation  coef¬ 
ficients  are  negative  for  extensive  ranges  of  the  separation  distance 
r^.  Under  these  conditions  interpretation  of  the  conventionally  defined 
scales  as  lengths  characteristic  of  the  average  size  of  the  turbulence^ 
eddies  is  open  to  question,  and  consideration  of  the  correlation  coeffi¬ 
cient  curves  is  advisable. 


The  correlation  coefficients  at  station  A  for  isotropic  turbulence 
with  the  spectrum  amplitude  function  G(k)  given  by  equation  (3l)  are 
obtained  from  equation  (7)  as 


(40) 


The  contour  integrations  used  to  obtain  equations  (40)  and  (4l)  are  not 
valid  when  r^  =  0;  hence  the  microscales  Ay  are  evaluated  from  the 

integral  relation  of  equation  (9) •  Such  evaluations  indicate  that 
AxA  «  X2A  -  o,  which  is  to  be  expected  in  view  of  the  neglect  of  vis¬ 
cosity  effects  in  the  analysis.  The  longitudinal  correlation  coefficient 
curve  of  equation  (40)  is  plotted  in  figure  8(a)  and  is  always  positive; 
the  lateral  correlation  coefficient  curve  of  equation  (41^  plotted  in 
figure  8(b)  reaches  its  zero  value  at  rxn  =2  (rx  =  2L]_  )  and  its 

minimum  value  at  r-^n  =  3  (r-g  =  3L-g  )  * 

The  downstream  correlation  coefficient  curves  (at  station  C)  have 
been  obtained  numerically  for  the  case  N  *  1  from  the  following 
rearrangement  of  equations  (7) : 
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In  evaluating  equations  (42)  and  (43),  values  of  the  integrand  were 
obtained  for  kp/n  ranging  from  0  to  50;  and  for  kp/n  greater  than 

Fr°(-)  G7) 

50,  in  view  of  the  asymptotic  behavior  of  the  functions  n 


the  integrand  was  approximated  as 


cos  kprp 
(kp/n)2 


f/(0)  (q^' 


Typical  downstream  longi¬ 


tudinal  and  lateral  correlation  coefficient  curves  (for  the  case 
Mg  =  0.05,  Me  =  2.00,  K  =  2,  N  =  l)  are  also  plotted  in  figures  8(a) 
and  8(b) ,  respectively,  to  indicate  the  changes  resulting  from  passage 
of  initially  isotropic  turbulence  through  a  given  screen  and  contraction. 
Although  the  downstream  lateral  correlation  coefficient  is  shown  in  fig¬ 
ure  8(b)  to  reach  slightly  negative  values,  it  is  believed  that  these 
are  the  result  of  unavoidable  "round-off"  errors  in  computation  of  the 
Fourier  transforms  and  that  the  coefficient  is  actually  always  positive, 
consistent  with  the  corresponding  spectrum  curve  of  figure  6(b),  which 
has  its  maximum  value  at  zero  wave  number. 


The  correlation  between  simultaneous  fluctuation  velocities  at  two 
points  a  distance  rp  apart  will  decrease  more  rapidly  with  increasing 
values  of  rp  when  the  eddies  comprising  the  turbulent  field  are  small 
than  when  the  eddies  are  large.  Figure  8(a)  thus  indicates  that  the 
longitudinal  scale  of  an  initially  isotropic  field  of  turbulence  is 
decreased  by  passage  through  the  particular  screen-contraction  config¬ 
uration  chosen.  Figure  8(b)  indicates  that  the  corresponding  lateral 
scale  is  increased. 

In  view  of  the  negative  values  attained  by  the  downstream  longi¬ 
tudinal  correlation  coefficient,  no  physical  meaning  can  be  assigned  to 
the  longitudinal  scale  ratio  defined  in  the  conventional  manner  by 
equation  (10).  For  example,  the  longitudinal  scale  ratio  reaches  a  zero 
value  even  though  the  longitudinal  turbulence  velocity  ratios  are  finite 
when  the  screen  pressure-drop  coefficient  NK  has  the  value  2.76.  The 
negative  values  attained  by  the  upstream  lateral  correlation  coefficient 
do  not  present  a  similar  anomaly  because  of  the  relation  between  the 
longitudinal  and  lateral  scales  in  the  case  of  isotropic  turbulence, 
namely,  LpA  =  2L2A. 
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The  difficulty  is  removed  if  an  effective  longitudinal  scale  Lq' 
is  defined  as  the  positive  area  under  the  corresponding  correlation 
curve.  Effective  longitudinal  scale  ratios  are  plotted  in  figure  9  and 
show  a  qualitative  similarity  with  the  conventional  ratios  shown  in  fig¬ 
ure  7(h).  For  a  given  value  of  the  screen  pressure-drop  coefficient 
NK,  the  effective  scale  ratio  decreases  with  increasing  speed  ratio  Zq 
to  a  minimum  value  at  Zq  =  27.4  for  which  the  contraction  parameter 

s?1  has  its  minimum  value.  For  a  given  contraction  the  effective  scale 
ratio  reaches  its  minimum  value  when  NK  Z  2.76. 


ESTIMATION  OF  DECAY  EFFECTS 

In  view  of  the  assumptions  of  inviscid  flow  and  small  fluctuation 
velocities  relative  to  the  main  stream,  the  preceding  analysis  is 
strictly  applicable  only  in  the  absence  of  the  turbulent  decay  processes 
(viscous  dissipation  and  turbulent  mixing) .  For  many  wind-tunnel  con¬ 
figurations,  effects  of  decay  upon  turbulence  are  of  the  same  order  of 
magnitude  as  the  screen-contraction  effects.  Correction  of  the  theo¬ 
retical  turbulence  velocity  ratios  may  therefore  prove  necessary  for 
practical  applications  of  the  theory. 

Selection  of  the  appropriate  decay  correction  presents  certain 
difficulties  inasmuch  as  there  is  a  lack  of  experimental  investigations 
of  axisyrnmetric  turbulence  decay.  Some  guidance  may  be  obtained  from 
the  theoretical  studies  of  Batchelor  (reference  4)  and  Chandrasekhar 
(reference  12)  on  axisyrnmetric  turbulence.  The  time  rates  of  change  of 
the  mean-square  velocity  components  are,  in  the  notation  of  reference  4: 

§t  1 )  “  *  4m0  +  2v("  10a  "  zlD  -  2c  -  14d) 

JL  =  2mQ  +  2v(-  10a  +  b  -  3d) 

In  these  equations  and  in  equations  (44)  and  (45),  the  symbol  V  rep¬ 
resents  the  kinematic  viscosity  coefficient.  The  corresponding  expres¬ 
sion  for  the  mean-square  resultant  velocity  is 

iL  (u-,2  +  2u92)  =  -  2t(30a  +  2c  +  20d)  (44) 

dt  x  1  c  ' 

For  isotropic  turbulence,  c  =  d  =  0  and  equation  (44)  becomes 

(“7  +  2“?)  =  (?“?)  -  -z',(30a) 


(45) 
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p  p 

The  velocity  components  u-^  and  u2  of  reference  4  are  identical  with 
q^2  and  q22  in  the  present  notation.  The  quantities  a,  b,  c,  and  d 

in  appropriate  groupings  represent  the  coefficients  in  the  series  expan¬ 
sions  in  r-^  for  the  longitudinal  and  lateral  velocity  correlation 
coefficients.  The  quantity  m^  depends  on  the  two-point  velocity- 
pressure  correlation  which  tends  to  zero  as  isotropy  is  approached.  For 
the  decay  of  isotropic  turbulence  in  a  constant-area  channel  during  the 
initial  period  wherein  both  inertia  and  viscous  forces  are  of  importance, 
equation  (45)  leads  to  the  semiempirical  relation  (reference  13) 


1 

3 


+  2 


0. 58t( Z-, ) 


(46) 


where  q^2  and  q^2  represent  the  mean-square  velocity  components  at 

any  station  downstream  of  the  reference  station  A  and  t (lj)  represents 
the  appropriate  decay  time. 


The  absence  of  the  velocity-pressure  correlation  term  niQ  in  both 
equations  (44)  and  (45)  suggests  that,  provided  the  quantity  (2c  +  20d) 
is  much  smaller  than  the  quantity  30a,  equation  (46)  may  yield  a  sat¬ 
isfactory  approximation  for  the  decay  of  the  mean-square  resultant  tur¬ 
bulent  velocity  in  axisymmetric  turbulence.  The  data  of  references  1 
and  14  tend  to  support  such  an  approximation.  The  assumption  that  the 
effects  of  the  screen-contraction  combination  and  the  decay  upon  the 
turbulent  velocity  ratios  proceed  independently  (see  reference  3)  leads 
to  the  relation 


Ul2)u  0  ^22)K 

j 

lql2)w  n  lq22)N 

Lc^)1  U2)aj 

3 

scd 

(?)AJ 

(47) 


where  the  subscript  sc  refers  to  the  turbulence  velocity  ratios 
obtained  in  the  absence  of  decay,  computed  from  equations  (25)  and  (26) 
and  listed  in  table  I,  and  the  subscript  scd  implies  that  the  effects 
of  initial  period  decay  have  been  included.  In  computing  J  from  equa¬ 
tion  (46),  the  decay  time  t(lj)  is  taken  as  the  time  required  for  a 
particle  at  local  main-stream  velocity  to  pass  through  the  screens  and 
the  contraction  starting  from  station  A.  This  implies  that  the  contrac¬ 
tion  affects  only  the  decay  time.  Some  question  exists  as  to  the  appli¬ 
cability  of  equation  (46)  and  hence  equation  (47)  for  damping  screens 
in  which  the  wire  diameters  are  usually  very  small. 
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A  comparison  of  the  theoretical  mean-square  resultant  turbulence 
velocity  ratio  corrected  for  decay  by  the  use  of  equation  (47)  with  the 
experimental  ratios  obtained  from  reference  1  is  shown  in  figure  10. 

The  mean-square  resultant  velocity  in  the  absence  of  decay  for  the  case 
of  a  single-screen-contraction  configuration  (W=l)  is  also  included  to 
show  the  magnitude  of  the  correction  involved  for  the  configuration  of 
reference  1.  The  .following  data  were  used  in  applying  the  decay  cor¬ 
rection:  UA  =  62.8  feet  per  second,  q^2 ) =  0.15  foot  per  sec¬ 
ond;  and  =  0.05  foot  (estimated).  The  screen  pressure-drop  coef¬ 

ficients  were  corrected  as  suggested  in  reference  14 


K  = 


K_ 


UA  dK 
~*°A 


where  Kg  designates  the  screen  pressure -drop  coefficient  measured  at  a 
given  speed  UA.  Although  the  single  experimental  points  obtained  for 
each  multiscreen  arrangement  do  not  check  the  decay  correction  as  well 
as  do  those  for  the  single-screen  arrangement,  the  limited  data  do  not 
warrant  any  refinement  of  the  correction  method  for  multiscreen- 
contraction  configurations. 

In  order  to  obtain  the  resolution  of  the  resultant  turbulence  veloc¬ 
ity  ratio  into  longitudinal  and  lateral  components,  some  knowledge  of 
the  velocity-pressure  correlation  is  required.  As  shown  in  reference  4 
the  effect  of  this  correlation  as  represented  by  the  term  jhq  is  to 
transfer  energy  from  the  larger  to  the  smaller  of  the  velocity  components, 
thus  providing  a  drive  towards  isotropy.  As  shown  in  table  I,  the  lon¬ 
gitudinal  component  will,  in  general,  be  much  smaller  than  the  lateral 
component  so  that  adjustment  of  the  longitudinal. component  is  more  crit¬ 
ical  than  adjustment  of  the  lateral  component.  The  magnitude  of  the 
longitudinal  component  is  governed  by  two  opposing  effects.  Turbulent 
decay  processes  reduce  this  component ;  the  drive  towards  isotropy  tends 
to  increase  it.  In  the  absence  of  any  quantitative  knowledge  concerning 
the  velocity-pressure  term  mg,  the  simplest  assumption  to  be  made  is 
that  the  longitudinal  turbulence  velocity  ratio  may  be  corrected , for 
decay  and  isotropy  drive  by  taking  an  average  of  the  values  for  zero 
decay  and  isotropic  decay  or 


Consistent  values  of  the  lateral  turbulence  velocity  ratio  are  then 
obtained  from  the  longitudinal  velocity  ratio  of  equation  (48)  and  the 
resultant  velocity  ratio  of  equation  (47). 
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The  comparison  shown  in  figure  11  provides  some  estimate  as  to  the 
agreement  that  might  "be  expected  between  the  predicted  turbulence- 
velocity-ratio  components  (corrected  for  decay)  and  the  experimental 
values.  The  agreement  shown  is  considered  satisfactory  for  most  engi¬ 
neering  applications.  The  theoretical  velocity  ratios  obtained  in  the 
absence  of  decay  are  included  for  the  case  N  =  1  to  indicate  the 
magnitude  of  the  correction. 

The  turbulence  scales  are  also  affected  by  the  turbulence  decay 
process,  tending  to  increase  as  the  decay  time  is  increased.  Under  the 
action  of  the  viscous  forces  the  smallest  eddies  are  dissipated  so  that 
the  average  eddy  size  (scale)  would  be  expected  to  increase.  For 
isotropic  turbulence,  the  change  in  scale  during  the  initial  period 
analogous  to  the  relation  given  for  the  fluctuation  velocity,  equa¬ 
tion  (46),  is 


Presumably  the  effect  of  decay  upon  the  scales  of  turbulence  could 
thus  be  obtained  by  a  procedure  similar  to  the  one  suggested  for  the 
fluctuation  velocities.  In  the  absence  of  any  experimental  data  such 
development  does  not  appear  warranted. 


CONCLUDING  REMARKS 

The  present  analysis  treats,  in  the  absence  of  turbulent  decay 
processes,  the  combined  effect  of  a  series  of  identical  damping  screens 
followed  by  a  stream  convergence  (or  divergence)  upon  the  mean-square 
fluctuation  velocities,  scales,  correlation  coefficients,  and  one¬ 
dimensional  spectra  of  a  field  of  turbulence  convected  by  a  main  stream. 
Numerical  results  are  presented  for  the  case  of  upstream  isotropic 
turbulence . 

The  limited  experimental  data  available  confirm  at  least  qualita¬ 
tively  some  of  the  theoretical  results  obtained  such  as  the  distortion 
of  an  initially  isotropic  field  of  turbulence  by  the  damping  screens  and 
stream  convergence  into  a  field  axisymmetric  about  the  main-stream 
direction  with  the  lateral  components  of  the  resultant  fluctuation  veloc¬ 
ity  larger  in  magnitude  than  the  longitudinal  component,  and  the  relative 
insensitivity  of  the  lateral  scale  of  turbulence  to  damping-screen  and 
stream-convergence  effects.  The  beneficial  effects  of  using  several 
screens  in  series  to  attain  a  given  over -all  screen  pressure-drop  coef¬ 
ficient  in  attenuating  the  fluctuation  velocities  are  also  substantiated. 
This  attenuation  is  accentuated  as  the  screen  coefficient  NK  is 
increased. 
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The  theory  predicts  certain  marked  changes  in  the  ordinates  of  the 
downstream  one -dimensional  spectra  and,  in  the  case  of  the  longitudinal 
spectra,  a  noticeable  distortion  of  shape  which  should  be  confirmable  by 
experiment.  The  longitudinal  downstream  correlation  coefficients  attain 
negative  values  over  a  large  range  of  the  separation  distance  r-^. 

Under  these  conditions,  the  scales  of  turbulence  as  conventionally 
defined  cannot  be  regarded  as  representative  of  the  average  eddy  size. 
Accordingly,  the  longitudinal  scales  have  been  redefined.  The  effect 
of  the  damping  screens  and  stream  convergence  is  to  decrease  the  lon¬ 
gitudinal  scale  and  to  increase  the  lateral  scale. 

An  approximate  method  of  correcting  the  predicted  turbulent  fluc¬ 
tuation  velocities  for  the  effects  of  turbulent  decay  is  presented. 
Tabulations  of  the  fluctuation  velocities  over  a  wide  range  of  conditions 
are  provided  for  convenience  in  engineering  applications. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  October  28,  1952 
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APPENDIX  A 

SYMBOLS 


The 

Al,2 

a2 

Bl,2,3 

b 

Cl,2,3 

Dl,2,3 

d 

El,2,3 


f 

G(k) 

g 

H 


following  symbols  are  used  in  this  report: 
parameter  groupings  defined  after  equation  (27) 
auxiliary  contraction  parameter,  a2  =  l/l  -  e 
parameter  groupings  defined  after  equation  (28) 

solidity  ratio  of  damping  screen 

parameter  groupings  defined  after  equation  (35) 

edge  lengths  of  volume  within  which  the  turbulence  field  is 
defined 

wire  diameter  of  damping  screen 

parameter  groupings  defined  after  equation  (36) 

Fp,  F2,  or  F5 

one -dimensional  longitudinal  spectral  density  (see  equation  (6)) 
one -dimensional  lateral  spectral  densities  (see  equation  (6)) 

auxiliary  wave-number  parameter,  f  3  s/rp  +  4/p^ 

amplitude  function  in  isotropic  spectrum  tensor  (see  equa¬ 
tions  (4)  and  (31)) 

auxiliary  wave-number  parameter,  g  s  s/|E  +  4a^/p^ 

constant  appearing  in  amplitude  function  of  special  isotropic 

spectrum  tensor,  H  =  ^  ( ll2/  (see  equation  (31)) 

■sr  n 


h 

i 


auxiliary  wave -number  parameter,  h 


y^r 
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J 

K 


k 


Lr 


^,0 


m 


N 


n 


P 

P 

Q 


turbulence  decay  factor  (see  equation  (46)) 

Ap 

screen  pressure-drop  coefficient,  K=  t  o 

=  PU^ 


amplitude  'of  vector  k:  k^  =  k-^  +  k2^  +  k3^ 


k  =  k..  ,  k„,  or  k,,;  wave-number  vector 
T  A.  La  O 


Li,  L2 ,  or  L3 

Lj  longitudinal  scale  of  turbulence  (see  equation  (10)) 

'  effective  longitudinal  scale  of  turbulence 

L2,  L3  lateral  scales  of  turbulence  (see  equation  (10)) 

stream  velocity  at  station  C  divided  by  stream  velocity  at 
station  B  (see  equation  (19)) 

stream  breadth  at  station  C  divided  by  stream  breadth  at  sta¬ 
tion  B  (see  equation  (19)) 

stream  height  at  station  C  divided  by  stream  height  at  sta¬ 
tion  B 


stream  Mach  number  at  station  B,  C 

mesh  designation  of  damping  screen  (reciprocal  of  center-to- 
center  distance  between  neighboring  wires) 

number  of  screens  in  series  (cascade) 

constant  appearing  in  amplitude  function  of  special  isotropic 
spectral  tensor,  n  =  1 


(L1)i 


constant 
static  pressure 

=  Q^,  Q2,  or  wave -amplitude  vector 


30 


NACA  TN  2878 


£  ly  =  Q.i>  0.2  >  or  03 i  turbulence -velocity-fluctuation  vector 

R^(r-j_)  correlation  coefficient  (see  equation  (7)) 

R^8(r)  correlation  tensor,  Ry§(r)  =  q.y(x)<l5(x+r)  . 

* 

r  rT  =  rl>  v2>  or  r3»  separation  vector 

O  .  p 

s  wave-number  parameter,  s  s  k^/y  +  1 

t  time 

t(Z1)  decay  time 


main-stream  velocity 

longitudinal  component  of  combined  turbulent  velocity  fluctua¬ 
tions  and  potential-flow  velocities  induced  by  screen 

longitudinal  root -mean-square  turbulence  velocity  ratio  (used 
in  table  I) 

lateral  root-mean-square  turbulence  velocity  ratio  (used  in 
table  I) 

lateral  components  of  combined  turbulent  velocity  fluctuations 
and  potential-flow  velocities  induced  by  screen 


Xy.  =  x-j_,  x2,  or  x3j  position  vector 


screen  deflection  parameter,  a  =  lim 

♦i-*o 


tan  \lfg 
tan  xlr-g 


,  ko  +  ki 


Fyg(k)  three-dimensional  spectral  tensor 


.  2  2  _  2  ,  2  _ 
4a,  cos  9  +  V  sin  9 

4  cos2  9  +  (i2  sin2  6 


Kronecker  delta;  5^.g  =  1  for  y  =  6  and  S^g  =  0  for  y  £  6 
axisymmetric  contraction  parameter,  g  = Ig2/^!2 

k22  +  k32 

2  -  ^ 

auxiliary  screen  parameter,  T]  =  — ^ - 

[xc  -  4 

polar  angle  (see  appendix  B) 

amplitude  of  vector  K_;  K2  =  Kq2  +  Kg2  +  Kg2 

K  =  ,  K„,  or  K  •  wave -number  vector  at  station  C 

y  1  2  3 

4a2k  2+V2^2 

- s - 5-—  (see  equation  (22a)) 

4k12  +  p2£2 

Aq,  A2,  or  A  3 

longitudinal  microscale  of  turbulence  (see  equation  (9)) 
lateral  microscales  of  turbulence  (see  equation  (9)) 


auxiliary  screen  parameter,  p  =  1  +  a  +  K 
auxiliary  screen  parameter,  v  =  1  +  a  -  aK 


auxiliary  screen  parameter, 


stream  density 


( u2 -a2 p2) kq2 

- 5 - 5—5-  (see  equation  (22b)) 

4kq2  +  P2C2 


main-stream  density  at  station  C  divided  by  main-stream  density 
at  station  B 


volume 
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cp 

X 

*1 
i2 

fl 

0) 

Superscripts : 

A  station  upstream  of  screens 

B  station  downstream  of  screens  and  upstream  of  contraction 

C  station  downstream  of  contraction 

*  complex  conjugate 

4 

Subscripts: 

A  station  upstream  of  screens 

<  K 

B  station  downstream  of  screens  and  upstream  of  contraction 

C  station  downstream  of  contraction 

N  number  of  like  ‘screens  in  series 

sc  only  effects  of  screens  and  contraction  present 

scd  effects  of  screen  and  contraction  corrected  for  initial  period 

of  decay 

1  longitudinal  component 

2,3  lateral  components 

/  . 


azimuth  angle  (see  appendix  B) 

Xy  =Xp>  "X»2>  or  ^3>  position  vector  (see  equation  (12)) 
angle  to  screen  normal  of  flow  incidence  upstream  of  screen 
angle  to  screen  normal  of  flow  emergence  downstream  of  screen 

k^d-c)2  £2  -  2k12(l-e)(ek12+£2) 

- - 2 — 2~2 -  (see  equation  (22d)) 

( ek-^  +£  ) 

ay  =  a>2>  or  053}  vorticity  vector 
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APPENDIX  B 


TURBULENCE  VELOCITY  AND  SCALE  RATIOS 
Velocity  ratios.  -  Using  spherical  polar  coordinates 

=  k  cos  0 

kg  =  k  sin  9  cos  cp 
k^  =  k  sin  9  sin  cp 
k2  —  kj2  +  kg2  +  kg2 

equations  (23)  and  (24)  may  he  put  in  the  form 


[rnc  Wl- 


Z22k2G(k) 


A^  sin2  0 


^1  (e  cos2  0  +  sin2  0)2 


[p22C(K)  +  r33C(jsON  =  a2[l22C(j5)  +  r33C(iii]N.1  + 


Z]G(k)  (yg-a2n2)  k2  A1^1  £2  sin2  9  cos2  8 
(4  cos2  0  +  n2  sin2  g)(e  cos2  8  +  sin2  8 )'c 


,  A  _  4a2  cos2  0  +  V2  sin2  0  / 

where  A  =  - - 5 - 5 - 5 - 

4  cos'2  0  +  \ic  sin^  0 

The  downstream  mean-square  fluctuation  velocities  are  given  hy 


^qT"2)  dKidK2dK3 


analogous  to  equation  (5).  Inasmuch  as  the  function  G(k)  appears  in 
the  expressions  for  the  energy  spectral  densities,  the  variable  of  inte¬ 
gration  will  be^changed  from  K,  to  k  so  that 

q^2  =  — ;L_.  jj  lr y-y-(K)  dkqdkgdkg.  Noting  that  .. 


ZqV  -co 


dkj_dkgdkg  =  k2sin  0  d0  dcp  dk,  the  downstream  mean-square  velocity  compo¬ 
nents  of  the  turbulent  field  are  obtained  from  equations  (Bl)  and  (B2)  as 
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S°kia(k)  “ 


an  sin5  e  ae 
(e  cos2  6  +  sin2 


and,  inasmuch  as  the  downstream  turbulence  will  be  axi symmetric  when 
the  upstream  turbulence  is  isotropic, 


WX  ■  M)n  ■  aZ  ML  +  e2(vz2-f^-1  jfW)  ak 


_ AN-1  sin3  6  cos2  0  d6 _ 

(4  cos2  0+^2  sin2  0)(c  cos2  0  +  sin2  0)2 


The  mean-square  velocity  components  of  the  upstream  isotropic  turbulence 
are  obtained  by  using  equation  (4)  as 


(q.12)  =  (q.22)  =  (l32)  =  Jo  k4G(k)  dk  '  dp  X  sin3@  d0  =  ^  Jo  k4G(k) 


dk 


The  turbulence  velocity  ratio  or  ratio  of  mean-square  fluctuation 
velocities  downstream  of  a  series  of  N  identical  screens  followed  by 
an  axisymmetric  contraction  to  the  corresponding  upstream  fluctuation 
velocities  is  then  given  for  the  longitudinal  and  lateral  components, 
respectively,  by 


/  AN  sin5  9  d9 

I  (a2  -  cos2  0)2 

n 


(B3) 


c  i 
•n  _ 

c 

h-1  ,  3(a2-l)2(v2-a2n2) 

j  '' n  A11*1  sin5  0  cos2  0  d0 

1 

WT  J 

Q  (4  cos20+|i2  sin2  0)(a2-cos2  0)2 

(B4) 


where 


a1- 


1 

1-  €* 


The  new  contraction  parameter 


is  introduced  for 


convenience  in  subsequent  calculations.  It  may  be  noted  that  the  veloc¬ 
ity  ratios  are  independent  of  the  amplitude  function  G(k)  which  appears 
in  both  the  isotropic  and  axisymmetric  spectral  tensors.  Equations  (B3) 
and  (B4)  appear  in  the  text  as  equations  (25)  and  (26),  respectively. 
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Turbulence  scale  ratios.  -  The  turbulence  scales  may  be  obtained 
from  the  energy  spectral  densities  as  indicated  by  equations  (6)  and 
(10).  Compatible  with  the  formulation 

\ 

CD 

(?)C  =  T^z'tt'frrrC{*-) 

—oo 


the  longitudinal  scale  at  station  C  is 


The  ratio  of  longitudinal  scale  downstream  of  a  series  of  N  identical 
screens  followed  by  an  axisymmetric  contraction  to  the  corresponding 
upstream  scale  or  longitudinal  scale  ratio  is  thus 
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The  corresponding  ratio  for  the  lateral  scales  is  obtained  in  a  similar 
manner  as 


These  relations  for  the  scale  ratios  do  not  require  that  the  upstream 
turbulence  be  isotropic.  Equations  (B5)  and  (B6)  appear  in  the  text  as 
equations  (37)  and  (38),  respectively. 


* 
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(f?C)  =  a?(jPC)  +  ^H(a2-l)2(v2-a2n2) 
\Z/N  N2'!!-!  l2  2 


‘  (toV+^eT  (ki2+^)V^ 
(4k12^??2)N(l£l2+^+?2)3[a2(kl2^2)  -  kl2J 


(C4) 


For  the  case  of  a  single-screen-axisymme trie -contract ion  configuration, 
integration  of  equations  (C3)  and  (C4)  yields  equations  (35)  and  (36)  of 
the  text,  respectively.  For  the  case  N  =  1,  the  quantity  (F2C)n  ^ 

designates  the  one-dimensional  lateral  spectrum  downstream  of  an  axisym- 
metric  contraction  in  the  absence  of  damping  screens  and  may  be  obtained 
from  reference  3. 
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TABLE  I.  -  TURBULENCE  VELOCITY  AND  SCALE  RATIOS 


ME ME , . 


'  m 


N 

K 

NK 

n 

H 

ID 

i(Vj£.2V22) 

Vl 

Mn 

N1 

NS 

N* 

1 
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0.0  0 

■an 

1.1  9  8  1  0  0 
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2.7  7  2  0  0  0 

0.4  8  36 

1 
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.0  2  3 

■Kftl 
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Continued.  TURBULENCE  VELOCITY  AND  SCALE  RATIOS 
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.331505 

.0  26  8 

.70  4  9 

.008098 

1.33  3  2 

l 

4.5  0 

4.5  0 

3.0  0  0 

3  5.8  6  6 

.000  64  7 

•45  1  72  4 

.30  136  5 

.0  25  4 

.672  1 

.01  3877 

1.3  3  31 

i 

6.0  0 

6.0  0 

3.0  0  0 

3  5.8  6  6 

.0005  01 

.3  54  92  6 

.23  6784 

.0  2  24 

.5  95  8 

.032835 

1.3  3  3  1 

1 

1  0.0  0 

1  0.0  0 

3.0  0  0 

3  5.8  6  6 

.00  0  32  5 

225862 

.15  0  6  8  3 

.0180 

.4  75  2 

.07  342  3 

1.3  3  32 

i 

0.0  0 

0  0.0  0 

4.0  0  0 

3  9.0  4  6 

0.0  0  3  8  4  4 

0.8  1  1  4 

0.5  4  22  0  0 

0.0  6  2  0 

0.9  0  0  6 

0.17  060  0 

1.3  3  30 

1 

.2  0 

0.2  0 

4.0  0  0 

3  9.0  4  6 

.003  2  30 

.734  08  7 

.490468 

.0  5  68 

.856  8 

.13  6  0  81 

1.3  3  2  9 

1 

.4  0 

0.4  0 

4.0  0  0 

3  9.0  4  6 

.002  7  42 

.664  1  3  2 

.44  36  6  9 

.05  2  4 

.814  9 

.10  718  5 

1.3  3  2  9 

i 

.6  0 

0.6  0 

4.0  0  0 

3  9.0  4  6 

•002  34  6 

.600693 

.40124  4 

.0  4  84 

.7  75  0 

.0  8  34  62 

1.3  3  2  9 

l 

1.0  0 

1.0  0 

4.0  0  0 

3  9.0  4  6 

.001747 

.4  90  79  1 

.32  7  77  6 

.0  418 

.700  6 

.04  8648 

1.3  329 

1 

1.5  0 

1.5  0 

4.0  0  0 

3  9.0  4  6 

.001277 

.3  92  67  3 

.26  2208 

.0  35  7 

.626  6 

.02  13  85 

1.3  3  2  9 

1 

2.0  0 

2.0  0 

4.0  0  0 

3  9.0  4  6 

.0009  93 

32722.7 

.21  848  2 

.0  315 

.57  2  0 

.0  0  66  5  4 

1.3  3  2  9 

1 

2.5  0 

2.5  0 

4.0  0  0 

3  9.0  4  6 

.000  8  09 

.28048  0 

.18  7  2  5  6 

.0  28  4 

.52  9  6 

.000676 

1.3  3  2  9 

l 

2.7  6 

2.7  6 

4.0  0  0 

3  9.0  4  6 

.0007  37 

.26108  5 

.17  4  3  0  2 

.0  27  1 

.5110 

.000000 

1.3  3  29 

1 

3.0  0 

3.0  0 

4.0  0  0 

3  9.0  4  6 

.0006  81 

.2  45  42  0 

.16  3  8  4  0 

.0  26  1 

.4  95  4 

.00  04  65 

1.3  3  2  9 

1 

3.5  0 

3.5  0 

4.0  0  0 

3  9.0  4  6 

.0005  8  9 

.218151 

.14  56  3  0 

.0  24  3 

.4  67  1 

.00  38  82 

1.3  3  29 

1 

4.0  0 

4.0  0 

4.0  0  0 

3  9.0  4  6 

.00  0  5  20 

.19  6  3  3  6 

.13  106  4 

.0  228 

.44  3  1 

.00  94  66 

1.3  3  2  9 

1 

4.5  0 

4.5  0 

4.0  0  0 

3  9.0  4  6 

.0004  67 

.1  7  8  4  8  7 

.11914  7 

.0  216 

.42  2  5 

.01  62  43 

1.3  3  2  9 

1 

6.0  0 

6.0  0 

4.0  0  0 

3  9.0  4  6 

•000361 

.1  4  0  2  3  9 

.09361  3 

.0  19  0 

.37  4  5 

.0384  77 

1.3  3  2  9 

i 

1  0.0  0 

10.0  0 

4.0  0  0 

3  9.0  4  6 

•000235 

.0  8  9  2  4.  3 

.05  9574 

.0153 

.2  987 

.08  54  75 

1.33  2  9 

1 

0.0  0 

0  0.0  0 

5.0  0  0 

4  0.8  3  5 

0-0  0  3  1  7  1 

0.3  4  8  0 

0.2  3  3  1  0  0 

0.0  5  6  3 

<15  8  9  $ 

0.18  9  20  0 

1.3  320 

l 

.2  0 

0.2  0 

5.0  0  0 

4  0.8  3  5 

.002647 

.3  14  84  7 

.210780 

.0515 

.56  11 

.15  17  9  5 

1.3  32  5 

l 

.4  0 

0.4  0 

5.0  0  0 

4  0.8  3  5 

.002  2  31 

.284  84  3 

.19  0  6  3  9 

.0  47  2 

.533  7 

.12  04  37 

1.3  32  5 

l 

.6  0 

0.6  0 

5.0  0  0 

4  0.8  3  5 

.0018  9  6 

.25  7  63  4 

.17  2  3  8  8 

.0435 

.5  07  6 

.0  9  4422 

1.3  3  25 

i 

1.0  0 

1.0  0 

5.0  0  0 

4  0.8  3  5 

.001395 

.2  104  9  8 

.14  0  7  9  7 

.0  37  4 

.45  8  8 

.055706 

1.3  3  25 

l 

1.50 

1.5  0 

5.0  0  0 

4  0.8  3  5 

.0010  07 

.1  6  8  4  1  5 

.112  6  12 

.0  317 

.4104 

.024805 

1.3  3  2  5 

l 

2,0  0 

2.0  0 

5.0  0  0 

4  0.8  3  5 

•000775 

.1  4  0  3  4  5 

.093822 

.0  27  8 

.3746 

.007799 

1.3  3  2  5 

i 

2.5  0 

2.5  0 

5.0  0  0 

4  0.8  3  5 

•0006  26 

.1  2  0  2  9  6 

.08  0406 

.0  25  0 

.34  6  8 

.000799 

1.3  3  2  5 

i 

2.7  6 

2.7  6 

5.0  0  0 

4  0.6  3  5 

•000568 

.1119  7  7 

.07  4841 

.0  238 

.33  4  6 

.000000 

1.3  3  2  5 

i 

3.0  0 

3.0  0 

5.0  0  0 

4  0.8  3  5 

•000524 

.10  5  2  5  8 

.0703  4  7 

.0  229 

.32  4  4 

.00  05  53 

1.33  25 

l 

3.5  0 

3.5  0 

5.0  0  0 

4  0.8  3  5 

•000451 

.09  3  56  3 

.06252  6 

.0  212 

.30  5  9 

.00  4634 

1.3  3  2  5 

l 

4.0  0 

5.0  0  0 

4  0.8  3  5 

.000  3  97 

.084  20  6 

.056270 

.0  199 

.2  90  2 

.011331 

1.3  3  2  5 

i 

4.5  0 

4.5  0 

5.0  0  0 

4  0.8  3  5 

•000  3  56 

.076  551 

.051153 

.0189 

.2767 

.01  9476 

1.3  3  2  5 

i 

6.0  0 

6.0  0 

5.0  0  0 

4  0.8  3  5 

•000  27  5 

.0601  4  7 

.04  019  0 

.0166 

.24  5  2 

.0  4  6  1  4  8 

1.3  3  25 

l 

1  0.0  0 

1  0.0  0 

5.0  0  0 

4  0.8  3  5 

.00018  2 

.0  38  275 

.025577 

.0135 

.19  5  6 

.10  13  81 

1.3  3  2  5 

i 

0.0  0 

0  0.0  0 

6.0  0  0 

419  16 

0-0027  09 

0.1  6  3  7 

0.1  1  0  0  0  0 

0.0  5  2  1 

0.4  0  4  6 

0.2  1  00  00 

1.3  3  10 

i 

.2  0 

0.2  0 

6.0  0  0 

4  1.9  1  6 

.002  24  5 

.14  8  10  7 

.099486 

.0474 

.3848 

.16  98  82 

1.3  3  16 

i 

.4  0 

0.4  0 

6.0  0  0 

4  1.9  1  6 

.001877 

133992 

.089954 

.0433 

.366  0 

.13  58  97 

1.3  3  16 

l 

.6  0 

0.6  0 

6.0  0  0 

4  1.9  1  6 

.001582 

.12119  3 

.081323 

.0  398 

.34  8  1 

.10  7  3  71 

1.3  3  16 

1.0  0 

1.0  0 

6.0  0  0 

4  1.9  1  6 

.00114  8 

.099  0  1  9 

.066395 

.0  3  39 

.3147 

.06  42  35 

1.3  3  16 

i 

1.5  0 

1.5  0 

6.0  0  0 

4  1.9  1  6 

•000  8  16 

•07922  2 

.053087 

.0  28  6 

.2815 

.02  90  36 

1.3  3  1  6 

1 

2.0  0 

2.0  0 

6.0  0  0 

4  1.9  1  6 

•000621 

.066  01  8 

.0  4  4  2  1  9 

.0  2  49 

.25  6  9 

.00  92  43 

1.3  3  17 

l 

2.5  0 

2.5  0 

6.0  0  0 

4  1.9  1  6 

.00  04  96 

.056  58  6 

.03  786  9 

.0  22  3 

.2  37  9 

.00  09  56 

1.3  3  17 

1 

2.7  6 

2.7  6 

6.0  0  0 

4  1.9  1  6 

.0004  4  9 

•052  67  3 

.035265 

.0  212 

.22  95 

.000000 

1.3  3  17 

1 

3.0  0 

3.0  0 

6.0  0  0 

4  1.9  1  6 

.000413 

.049  51  3 

.03  3146 

.0  20  3 

.2225 

.000666 

1.3  3  17 

1 

3.5  0 

3.5  0 

6.0  0  0 

4  1.9  1  6 

.000354 

.044  011 

.02  9459 

.0188 

.20  9  8 

.005610 

1.3  3  1  7 

j 

4.0  0 

4.0  0 

6.0  0  0 

4  1.9  1  6 

.00  0  3  11 

•039  609 

.02651  0 

.0176 

.19  90 

.01  3756 

1.3  3  1  7 

1 

4.5  0 

4.5  0 

6.0  0  0 

4  1.9  1  6 

.000  27  8 

.036  00  8 

.02  4098 

.0  167 

.18  98 

.023677 

1.3  3  1  7 

l 

6.0  0 

6.0  0 

6.0  0  0 

4  1.9  1  6 

.0002  15 

.0  28  292 

.01  8933 

.0  14  7 

.16  8  2 

.05  6027 

1.3  3  1  7 

1 

1  0.0  0 

1  0.0  0 

6.0  0  0 

4  1.9  1  6 

.000144 

.0  18  00  3 

.012050 

.0120 

.13  4  2 

.12  1.0  0  8. 

1.3  3  1  8 

l 

0.0  0 

0  0.0  0 

7.0  0  0 

4  2.6  1 1 

0.002  363 

0.0  8  3  6  6 

0.0  5  6  60  0 

0.0  4  8  6 

0.2  8  9  2 

0.2  3  320  0 

1.3  3  1  6 

1 

.2  0 

0.2  0 

7.0  0  6 

4  2.6  1 1 

.0019  39 

.075  716 

.051124 

.0440 

.2  75  2 

.19  03  0  3 

1.3  3  0  1 

1 

.4  0 

0.4  0 

7.0  0  0 

4  2.6  1 1 

.001  6  07 

.068  4  99 

.04  6202 

.0  4  01 

.2617 

.15  3  603 

1.3  3  0  1 

1 

.6  0 

0.6  0 

7.0  0  0 

4  2.6  1 1 

.001  343 

.061  95  5 

.04  1751 

,0  366 

.24  89 

.12  2  4  0  2 

1.3  3  0  1 

1 

1.0  0 

1.0  0 

7.0  0  0 

4  2.6  1 1 

.00  09  6  0 

.05061  9 

.0  3  4  0  6  6 

.0310 

.2250 

.074378 

1.3  3  0  2 

1 

1.5  0 

1.5  0 

7.0  0  0 

4  2.6  1  1 

.00067  1 

.04  049  8 

.02  722  2 

.0  2  59 

.20  12 

.0  3  4200 

1.3  3  0  2 

1 

2.0  0 

2.0  0 

7.0  0  0 

4  2.6  1 1 

.000  5  03 

.033  74  8 

.0  2  2  6  6  6 

.0  224 

.183  7 

.011043 

1,3  3  0  2 

l 

2.5  0 

2.5  0 

7.0  0  0 

4  2.6  1  1 

.00  0  39  8 

.0  28  92  6 

.019417 

.0  199 

.170  1 

.001155 

1.3  3  0  3 

l 

2.7  6 

2.7  6 

7.0  0  0 

4  2.6  11 

.0003  58 

.0  26  92  5 

.018069 

.0189 

.16  41 

.000000 

1.3  30  3 

1 

3,0  0 

3.0  0 

7.0  0  0 

4  2.6  IX 

.00  0  32  8 

.0  25  3  0  9 

.0  1  6  9  8  2 

.0181 

.15  9  1 

.00  0811 

1.3  30  3 

1 

3.5  0 

3.5  0 

7.0  0  0 

4  2.6  1 1 

.000280 

.022  4  9  7 

.01  5091 

.0167 

.15  0  0 

.00  6861 

1.3  3  0  3 

1 

4.0  0 

4.0  0 

7.0  0  0 

4  2s6  1 1 

.000245 

.02024  7 

.01  3580 

.0157 

.14  2  3 

.01  68  69 

1.3  3  0  4 

1 

4.5  0 

4.5  0 

7.0  0  0 

4  2.6  1 1 

.000219 

.0  18  40  6 

.012344 

.0148 

.135  7 

.02  90  61 

1.3  3  0  4 

1 

6.0  0 

6.0  0 

7.0  0  0 

4  2.6  1 1 

.000170 

.0  14  461 

.009697 

.0130 

.12  0  3 

.06  84  95 

1.3  3  0  5 

1 

1  0.0  0 

10,0  0 

7.0  0  0 

4  2.6  1 1 

.0001  17 

.0  09  201 

.006173 

.0108 

.0  95  9 

.14  4  4  16 

13  30  6 

2662 
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TABLE  I.  -  Continued.  TURBULENCE  VELOCITY  AND  SCALE  RATIOS 


N 

K 

NK 

Mc 

vi2 

V22 

i(V12+2V22) 

vi 

V2 

MS 

"N* 

(Mfi 

MO1 

3 

3 

3 

3 

3 

3 

0.2  0 
.4  0 
.6  0 

1.0  0 

1.5  0 

2^0  0 

0.6  0 

1.2  0 

1.8  0 

3,0  0 

4.5  0 

6.0  0 

0.0  2  3 
.0  2  3 
.0  2  3 
.0  2  3 
.0  2  3 
.0  2  3 

0.4  6  4 
.4  6  4 
.4  6  4 
.4  6  4 
.4  6  4 
.4  6  4 

0.5  4  0  0  4  0 
.17  9  0  71 
.062658 
.010  8  2  9 
.003090 
.0016  02 

0.5  0  6  0  6  7 
.298  411 
.19  13  8  2 
.09  0  657 
.04  3429 
.0  24  431 
.0  10061 

0.5  1  7  3  9  1 
.258631 

4.4  847  4 
.064048 
.029982 
.016821 
.006926 

0.7  3  4  9 
.4  232 
.2  503 
4.041 
.0  5  56 
.0400 
.0  25  6 

0.7  1 1  4 
.54  6  3 
.4  37  5 
.30  11 
.208  4 
.15  6  3 
.1003 

3.5  8  5  6  8  7 
3.332176 

1.9  6  9  1 0  0 
.932930 
.10  8  4  6  2 
.002965 
.000001 

0.6  7  9  5 
.8  5  33 
.98  4  5 

1.1  3  3  5 

1.2  1 1  9 

1.3  4  6  7 

12  7  71 

_ 

3 

3 

3 

3 

3 

3 

0.2  0 
,4  0 
.6  0 

1.0  0 

1.5  0 

2.0  0 

0.6  0 

1.2  0 

1.8  0 

3.0  0 

4.5  0 

6.0  0 

0.0  4  0 
.0  4  0 
.0  4  0 
.0  4  0 
.0  4  0 
.0  4  0 

0.8  0  0 

.8  0  0 

.8  0  0 

.8  0  0 

.8  0  0 

.8  0  0 

0.4  1  4  5  5  9 
.15  0  9  15 
.0579  97 
.010863 
.002433 
.0009  4  8 

0.6  0  9  2  9  6 
.414814 
.29058  8 
.1  4  9  0  4  5 
.073  531 
.041  722 
017244 

0.5  4  4  3  8  4 
.326848 
.213058 
.10  2  9  8  4 
.049832 
.028130 

01  1595  ... 

0.6  4  3  9 
.3  885 
.2  4  08 
.104  2 
.0  4  93 
.0  308 

.0172 

0.7  8  0  6 
.64  41 
.5391 
.3861 
.27  12 
.20  4  3 
1313 

.931719 
.716250 
.3130  94 
.046378 
.00  1686 
on  00  01... 

l!0  5  7  9 
1.117  4 

1.1 8  8  2 

1.2  3  3  5 

1.2  5  8 1 

1.2  8  4  2 

-3.  _ 

3 

3 

3 

3 

3 

3 

0.2  0 
.4  0 
.6  0 

1.0  0 

1.5  0 

2.0  0 

0.6  0 

1.2  0 
1.0  0 
3.0  0 
4.5.0 
6.0  0 

0.0  6  0 
.0  6  0 
.0  6  0 
.0  6  0 
.0  6  0 
.0  6  0 

0  1.2  0  0 

1.2  0  0 

1.2  0  0 

1.2  0  0 

1.2  0  0 

1.2  0  0 

0.3  3  4  2  8  1 
.13  8  2  4  5 
.061227 
.014  9  01 
.003037 
.001410 
.000348 

0.774  595 
.5571  92 
.403  98  2 
.214  731 
.10  7  9  3  5 
.06174  8 
..025  71  5 

0.6  2  7  8  2  4 
.417543 
.289731 
4.4  8121 
.073236 
.041635 

■Q1725? 

0.5  7  8  2 
.3  718 
.2  47  4 
.12  21 
.0  619 
.0  37  5 
.0187 

0.8  8  0  1 
.74  6  5 
.6  35  6 
.46  3  4 
.32  8  5 
.24  8  5 
.160  4 

.452149 
.301609 
.10  1467 
.0130  76 
.000504 
.000000 

1.1  8  0  6 

1.2  0  4  8 

1.2  3  6  3 

1.2  5  9  7 

1.2  7  4  3 

1.3  9  0  9 

3 

3 

3 

3 

3 

3 

0.2  0 
.4  0 
.6  0 

1,0  0 

1.5  0 

2.0  0 

0.6  0 
1.2  0 
1.8  0 
3.0  0 
4.5  0 
6.0  0 

0.10  0 
.10  0 
.10  0 
.10  0 
.10  0 
.10  0 

0  2.0  0  0 

2.0  0  0 

2.0  0  0 

2.0  0  0 

2.0  0  0 

2.0  0  0 

0.2  33  8  20 
.114  5  3  5 
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Figure  1.  -  Configuration  treated  in  analysis. 
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Figure  2.  -  Action  of  damping  screen  on  components  of  combined  turbulent  and 

induced  velocities  at  screen. 
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Figure  4.  -  Variation  of  root-mean-square  turbulence  velocity  ratio  with  speed  ratio  (Mg  of  0.05)  and  screen  pressure-drop  coef¬ 
ficient  K  in  absence  of  turbulence  decay  for  single-screen-axisymmetrlc-contraction  configurations  with  upstream  isotropic 
turbulence. 
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Speed  ratio, 

(b)  Lateral  turbulence  velocity. 

Figure  4.  -  Concluded.  Variation  of  root-mean-square  turbulence  velocity  ratio  with  speed  ratio  (Mg  of  0.05)  and  screen  pressure 
drop  coefficient-  K  in  absence  of  turbulence  decay  for  single-screen-axisymmetric-contraction  configurations  with  upstream  iso¬ 
tropic  turbulence. 
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Figure  5.  -  Effect  of  multiple  screens  N  and  speed  ratio  (Mb  of  0.05)  on  root-mean-square  turbulence  velocity  ratio  in  absence 
of  turbulence  decay  for  screen-axisymmetric-contraction  configurations  with  upstream  isotropic  turbulence  and  constant  screen 
losses.  Over-all  screen  pressure-drop  coefficient,  NK,  6. 


absence  of  turbulence  decay  for  screen-axlsymine trie -contraction  configurations  with  upstream  Isotropic  turbulence  and  constant 
screen  losses.  Over-all  screen  pressure-drop  coefficient,  NK,  6. 
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(a)  Longitudinal  spectra. 

Figure  6.  -  Comparison  of  one-dimensional  spectra  in  absence  of  turbulence  decay  fo 
and  for  single-screen-contraction  configurations  for  upstream  isotropic *turbulenc 
amplitude  function  G(k)  =  H(k2  +  n2)-3.  MB,  0.05;  Mc,  2.00;  llf  29.822. 
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(b)  Longitudinal  scale  ratio.  Longitudinal  scale  ratio  equals  zero  for  K*  2.76. 

Figure  7.  -  Concluded.  Variation  of  scale  ratio  with  speed  ratio  (Mg  of  0.05)  and  screen  pressure- 

drop  coefficient  K  in  absence  of  turbulence  decay  for  s ingle -screen-axi symmetric- contract ion 
configurations  with  upstream  isotropic  turbulence. 
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(a)  Longitudinal  correlation  coefficients. 

Figure  8.  -  Comparison  of  correlation  coefficients  in  absence  of  decay  for  a 
screen-contraction  configuration  (Mg  =  0.05,  Mq  =  2.00,  K  =  2,  N  =  l)  with 

upstream  isotropic  turbulence  having  amplitude  function  G(k)  =  H(k2  +  n2 
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Figure  8.  -  Concluded.  Comparison  of  correlation  coefficients  in  absence  of 
decay  for  a  screen- contraction  configuration  (Mg  =  0.05,  Mq  =  2.00,  K  •=  2, 

N  =  l)  with  upstream  isotropic  turbulence  having  amplitude  function 
G(k)  =  H(k2  +  n2)-3. 
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Speed  ratio, 

Figure  9.  -  Variation  of  effective  longitudinal  scale  ratio  in  absence  of  turbulence  decay  with  speed  ratio  (Mg  of  0.05)  and  screen 
pressure-drop  coefficient  K  for  single-screen-axlsymmetric-contraction  configuration  with  upstream  isotropic  turbulence  having 
amplitude  function  G(k)  =  H(k2  +  n2)'3. 


Figure  10.  -  Comparison  of  theoretical  mean-square  resultant  turbulence  velocity  ratios  corrected  for  decay  with  experiment  of 
reference  1.  Speed  ratio  6.7  (Mg  =  0.05,  M^i  =  0.34);  N  screens  in  series;  upstream  isotropic  turbulence;  scale 

.1^ )  0.05  foot  (estimated). 
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Figure  11.  -  Comparison  of  theoretical  root -mean- square  longitudinal  and  lateral  turbulence  velocity  ratios  corrected  for  decay  with 
eiqjeriment  of  reference  1.  Speed  ratio  l±,  6.7  (Mg  =  0.05,  Mq  =  0.34);  N  screens  in  series;  upstream  isotropic  turbulence; 
scale  LgA,  0.05  foot  (estimated).  ✓ 


